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Abstract 

The DGLAP analyses of the nuclear partem distribution functions (nPDF) based 
on the global fits to the data are reviewed, and the results from EKS98 and 
HKM are compared. The usefulness of measuring hard probes in pA colli- 
sions, at the LHC in particular, is demonstrated. 

1. Introduction 

Inclusive cross sections for hard processes A+B — > c+ X involving a sufficiently large scale Q ^> Aqcd 
are computable by using collinear factorization. In the leading twist approximation power corrections 
~ l/Q 2 are neglected and 



da(Q 2 , \/~s)ab^c+x = 

i,j=q,q,9 



Z A ff /A ( Xl ,Q 2 ) + (A- Z A )f^(x u Q 



n/A, 







Z B ^ /B (X2, Q 2 ) + {B- Z B )f" /B (x 2 , Q 2 : 



da(Q 2 ,xi,x 2 ] 



ij—>c-{-x 



(1) 



where A, B are the colliding hadrons or nuclei containing Z A and Zg protons correspondingly, c the 
produced partem, x and X anything, da{Q 2 ,x\,X2)ij^ c +x the perturbatively calculable differential cross 
section for the production of c at the scale Q, and xi 2 ~ Q/y/s the fractional momenta of the colliding 
partons i and j. The number distribution function of the parton flavour i of the protons (neutrons) in A 
is denoted as /f {f^ A )> an d similarly for partons j in B. 

In the leading twist approximation multiple scattering of the bound nucleons does occur but all 

collisions are independent, correlations between partons from the same object A are neglected, and only 

j>/A 

one-parton densities are needed. The parton distribution functions (PDF) f[ are universal quantities 
and applicable in all collinearly factorizable processes. The PDF cannot be computed by perturbative 
methods, and so far it has not been possible to compute them from first principles, either. Thus, nonper- 
turbative input from data on various hard processes is needed for the extraction of the PDF. However, 
once the PDF are known at some initial (lowest) scale Qo ^> Aqcd> the QCD perturbation theory 
predicts the scale evolution of the PDF to other (higher) values of Q 2 in form of the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equations [ 1 1. 

The method to extract the PDF from experimental data is well established in the case of the free 
proton: the initial (non-perturbative) distributions are parametrized at some Qq, and evolved to higher 
scales according to the DGLAP equations. Comparison with the data is made at various regions of 
the (x, Q 2 ) -plane, and the parameters of the initial distributions ff (x, Qq) become fixed when the best 
global fit is found. The data from deeply inelastic Ip scattering (DIS) are of main importance in these 
global DGLAP fits, especially the HERA data at small values of x and Q 2 . The sum rules for momentum, 
charge and baryon number give further constraints. In this way, through the global DGLAP fits, groups 
like MRST Q, CTEQ Q or GRV Q| obtain their well-known sets of the PDF of the free proton. 
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The nuclear parton distribution functions (nPDF) differ in magnitude from the PDF of the free 
proton. In the measurements of the structure function F 2 = Z&F 2 A + (A — Za)F 2 ' A of nuclear 
targets in DIS (see e.g. [5] for references) and especially of the ratio 

A = \da lA /dQ 2 dx ^ ^F A (x,Q 2 ) 

H F2 (X, ^ ) - l d(jlD/dQ2dx ~ l p D^ Q2) , W 

the following nuclear effects have been discovered as a function of Bjorken-x: 

• shadowing; a depletion Rp < 1 at x < 0.1, 

• anti-shadowing; an excess Rp 2 > 1 at 0.1 < x < 0.3, 

• EMC effect; a depletion at 0.3 < x < 0.7, and 

• Fermi motion; an excess towards x — ► 1 and beyond. 

The Q 2 dependence of Rp 2 is weaker and has thus been more difficult to measure. Data with high 
enough precision, however, exist: NMC has some years ago discovered a clear Q 2 dependence in the 
ratio da^/da^ [6], i.e. the scale dependence of the ratio F 2 Sn /F 2 c , at x > 0.01. Since F\ 



V(n)/A 
2 



J2 q e q x [fq + /f ] + 0(a s ), the nuclear effects in the ratio Rp 2 directly translate into nuclear 
effects in the parton distributions; /f ^ ff. 

The nPDF /f also obey the DGLAP equations in the large-Q 2 limit, and they can be determined 
by using a similar global DGLAP fit procedure as in the case of the PDF of the free proton. Pioneering 
studies of the DGLAP evolution of the nPDF are e.g. 0|HJ|9jnP |. References for various other studies of 
perturbative evolution of the nPDF and also to simpler Q 2 -independent parametrizations of the nuclear 
effects in the PDF can be found e.g. in lHT1IT2l . The nuclear case is, however, more complicated because 
of additional variables, the mass number A and the charge Z, and, because the number of data points 
available in the perturbative region is more limited than for the PDF of the free proton. The DIS data 
plays the dominant role in the nuclear case as well. However, as illustrated by Fig.Q no data are available 
from nuclear DIS experiments below x < 5 • 10~ 3 at Q 2 > 1 GeV 2 . This makes the determination of the 
nuclear gluon distributions especially difficult. Further constraints for the global DGLAP fits of the nPDF 
can be obtained from e.g. the Drell-Yan (DY) process measured in fixed-target pA collisions JT3lfT4l . 
Currently, there are two sets of nPDF available which are based on the global DGLAP fits to the data: (i) 
EKS98 CUED (the code in HUGH), and (ii) HKM (19) (the code in l20l V We shall compare the main 
features of these two analyses and comment on their differences below. 

2. Comparison of EKS98 and HKM 

2.1 EKS98 and overview of constraints available from data 

The parametrization EKS98 prepared in [15] is based on the results of the DGLAP analysis in |16J 
and its follow-up in fi"5ll . We shall here refer to these together as EKRS. In the EKRS approach the 
nPDF, the parton distributions of the bound pn 
corresponding distributions in the free proton, 



v/A 

nPDF, the parton distributions of the bound protons, n' , are defined through the modifications of the 



RA( X Q2\ = fi_ ( X >Q 2 ) (3) 

where the PDF of the free proton are from a known set, such as CTEQ, MRS or GRV. As in the case 
of the free nucleons, for isoscalar nuclei the parton distributions of bound neutrons are obtained through 
isospin symmetry, = f^!t and f™!t = - Although exact only for isoscalar and mirror nuclei, 
this is expected to be a good first approximation for all nuclei. 

To simplify the determination of the input nuclear effects for valence and sea quarks, the following 
flavour-independent initial modifications are assumed: R a v (x,Qq) R A y (x, Qq) ^ Ry{x,Q 2 ) ), and 




X 



Fig. 1: The average values of x and Q 2 of the DIS data from NMC EIISIE) (triangles) and E665 EHUD (diamonds) in 
I A, and of X2 and M 2 of the DY dilepton data 1 13 1 (squares) in pA. The heavy quark mass scales are shown by the horizontal 
dashed lines. The initial scale Qq is m 2 in EKRS and 1 GeV 2 in HKM. The rest of the figure: see the text in Sec. 13. II 

R£(x, Qq) sa Rj(x, Qq) ~ Rf(x, Qq) pa Rs(x, Qq). Thus only three independent initial ratios, Ry, 
Rg and Rq are to be determined at Qq = m 2 = 2.25 GeV 2 . Note also that the approximations above are 
needed and used only at Qq in 1161 . In the EKS98-parametrization [ 15 1 of the DGLAP-evolved results, 
it was observed that to a good approximation R^ v = R^ and Rf = Rf for all Q 2 . Further details of 
the EKRS analysis can be found in [ 16 1, here we summarize the constraints available in each region of 
x. Consider first quarks and antiquaries: 

• At x > 0.3 the DIS data constrains only the ratio Ry : valence quarks dominate F^, so Rp sa R v 
but Rg and Rq are left practically unconstrained. An EMC effect is, however, assumed also for 
the initial Rg and Rq since in the scale evolution the EMC effect of valence quarks is transmitted 
to gluons and further on to sea quarks (9). In this way the nuclear modifications Rf also remain 
stable against the evolution. 

• At 0.04 < x < 0.3 both the DIS and DY data constrain Rg and R v but from different regions of 
Q 2 , as shown in Fig. ^ In addition, conservation of baryon number restricts Ry. The use of DY 
data [ 13 ] is essential in order to fix the relative magnitude of R v and Rg, since the DIS data alone 
cannot distinguish between them. As a result, no antishadowing appears in Rg(x, Qq). 

• At 5 • 1CT 3 < x < 0.04 only DIS data exist in the region Q > 1 GeV where the DGLAP analysis 
can be expected to apply. Once Ry is fixed by the DIS and DY data at larger x, the magnitude of 
nuclear valence quark shadowing in the EKRS approach follows from the conservation of baryon 
number. As a result of these contraints, nuclear valence quarks in EKRS are less shadowed than 
the sea quarks, Ry > Rg. 

• At x < 5 • 10~ 3 , as shown by Fig. ^ the DIS data for the ratio Rp lie in the region Q < 1 GeV 
where the DGLAP equations are unlikely to be applicable in the nuclear case. Indirect constraints, 
however, can be obtained by noticing: (i) At the smallest values of x (where Q <C 1 GeV), Rp 2 
depends only very weakly, if at all, on x 1221 1241 . (ii) 5(F 2 Sn /F 2 c )/aiog Q 2 > at x ~ 0.01 0, 
indicating that 3R f a / d log Q 2 > 0. (iii) Negative log Q 2 -slopes of Rp 2 at x < 5 • 10 -3 have not 



been observed [22|. Based on these experimental facts, it is assumed in the EKRS approach that 
the sign of the log<2 2 -slope of Rp 2 remains non-negative at x < 0.01 and therefore the DIS data 
1221 1241 in the non-perturbative region gives a lower bound for Rp 2 at Qq at very small x. The sea 
quarks dominate over the valence quarks at small values of x, so only Rg becomes constrained by 
the DIS data; R$ remains restricted by baryon number conservation. 

Pinning down the nuclear gluon distributions is difficult in the absense of stringent direct constraints. 
What is available in different regions of x can be summarized as follows: 

• At x > 0.2 no experimental constraints are currently available for the gluons. Conservation of 
momentum is used as an overall constraint but this alone is not sufficient to determine whether an 
EMC effect exists for gluons or not: only about 30 % of the gluon momentum comes from x > 0.2, 
so fairly sizable changes in R G (x, Qq) in this region can be compensated by smaller changes at 
x < 0.2 without violating the constraints discussed below. As mentioned above, in the EKRS 
approach an EMC effect is initially assumed for R G . This guarantees a stable scale evolution; the 
EMC effect remains there also at larger Q 2 . 

• At 0.02 < x < 0.2 the Q 2 dependence of the ratio F 2 Sn /F 2 c measured by NMC |6) sets the 
currently most important constraint for R G , as first pointed out in l26l . In the small-x region 
where gluons dominate the DGLAP evolution, the Q 2 dependence of F 2 (x, Q 2 ) is dictated by 
the gluon distribution as H3 dF 2 (x,Q 2 )/dlogQ 2 ~ a s xg(2x,Q 2 ). This leads to fT6l fTTl 
dR F2 (x,Q 2 )/dlogQ 2 ~ a s [R G (2x,Q 2 ) - Rp 2 (x,Q 2 )]xg(2x,Q 2 )/F 2 D (x,Q 2 ). The log Q 2 
slopes of F^ n /F^ measured by NMC [6] therefore constrain Rq. Especially, as discussed in detail 
in EDEH the positive log Q 2 -slope of F 2 Sn /F 2 c measured by NMC indicates that R G {2x, Q§) > 
Rp 2 (x, Qq) at x ~ 0.01. Thus, within the DGLAP framework, a gluon shadowing at x ~ 0.01 
which would be much stronger than the shadowing of antiquarks, such as suggested e.g. in [28 1, 
is not supported by the NMC data. The antishadowing in the EKRS gluons follows from the con- 
straint R G (0.03, Qq) « 1 imposed by the NMC data (see also [26 1) combined with the requirement 
of momentum conservation. The EKRS antishadowing is consistent with the E789 data [29| on 
D-meson production in pA collisions (notice, however, the large error bar of the data point), and 
it seems to be supported by the J/ty production in DIS, measured by NMC 1 30 1 . 

• At x < 0.02, stringent experimental constraints do not exist for the nuclear gluons at the moment. 
It should be emphasized, however, that the initial R G in this region is directly connected with the 
initial Rp 2 . As discussed above, related to quarks at small x, taking the DIS data on Rp 2 in the 
non-perturbative region as a lower limit for Rp 2 at Qq corresponds to 8R f a / d log Q 2 > 0, and 
thus R G {x < 0.02, Q 2 q) > Rp 2 (x/2,Ql) > Rp 2 (x/2,Q 2 < lGeV 2 ). The observation in ED 
was that setting R G (x, Qq) Rp 2 (x, Qq) at x < 0.01 fulfills this constraint. This approximation 
remains fairly good even after the DGLAP evolution from Qq ~ 1 GeV to Q ~ 100 GeV, see [ 15 1. 

As explained in detail in [ 16 1, in the EKRS approach the initial ratios Ry(x, Qq), Rg(x, Qq) and 
R G (x, Qq) are constructed piecewize in different regions of x. Initial nPDF are computed at Qq, LO 
DGLAP evolution to higher scales is performed and comparison with DIS and DY data is made. The 
parameters in the input ratios are iteratively changed until the best global fit to the data is achieved. The 
determination of the input parameters in EKRS has so far been done only manually, the best overall fit is 
determined by eye. For the quality of the obtained fit, see the detailed comparison with the data in Figs. 
4-10 of [ 16 1. The parametrization EKS98 [ 15 1 of the nuclear modifications Rf (x, Q 2 ) was prepared on 
the basis of the results in 1 16 1. It was also shown that when the PDF of the free proton where changed 
from GRVLO [31 1 to CTEQ4L [32| (differing from each other considerably), the changes induced to 
Rf(x, Q 2 ) were within a few percents fi"5l . Therefore, accepting this range of uncertainty, the EKS98 
parametrization can be used together with any (LO) set of the PDF of the free proton. 



2.2 The HKM analysis 

In principle, the definition of the nPDF in the HKM analysis [19| differs slightly from that in EKRS: 
instead of the PDF of the bound protons, HKM define the nPDF as the average distributions of each 



flavour i in a nucleus A: f A (x, Q 2 ) = {Z/A)ff' (x, Q 2 ) + (1 - Z/A)fi 1 (x, Q 2 ). Correspondingly, 



EKS98 modifications represent average modifications. Flavour-symmetric sea quark distributions are 
assumed in HKM, whereas the flavour asymmetry of the sea quarks in EKRS follows from that of the 
free proton. 

An improvement relative to EKRS is that the HKM method to extract the initial modifications 
Wi(x, A, Z) at Qq is more automatic and more quantitative in the statistical analysis; the HKM analysis 
is strictly a minimum-x 2 fit. Also, with certain assumptions of a suitable form (see [ 19 1) for the initial 
modifications W{, the number of parameters has been brought conveniently down to seven. The form 
used is 



where Hi(x) = biX + qx 2 (analysis with a cubic polynomial is also performed, with similar results). 
Due to the flavour-symmetry, the sea quark parameters are identical for all flavours. For valence quarks, 
conservation of charge Z (not used in EKRS) and baryon number A are required; this fixes a Uv and a<j v . 
In the case of non-isoscalar nuclei w Uv ^ Wd v - Also momentum conservation is imposed; this fixes a g . 
Taking 0y, by, cy to be the same for uy and dy, and Pq = P g = 1, b g = —2c g , the remaining seven 
parameters by, cy, (3y, aq bq, cq, c g are determined by a global DGLAP fit which minimizes the x 2 - 

2.3 The comparison 

Irrespective of whether the best fit is found automatically or by eye, the basic procedure to determine 
the nPDF in the EKRS and HKM analyses is the same. The results obtained for the nuclear effects are, 
however, quite different, as can be seen in the comparison at Q 2 = 2.25 GeV 2 shown in Fig. |2] The 
main reason for this is that different data sets are used: 

• The HKM analysis [ 19 1 uses only the DIS data, whereas EKRS include also the DY data from pA 
collisions. As explained above, the DY data is very important in the EKRS analysis in fixing the 
relative modifications of valence and sea quarks at intermediate x. Preliminary results reported 
in ll33l show that when the DY data is included in the HKM analysis, the antiquark modifications 
will become more similar to those in EKRS. 

• The NMC data set on F A / F 2 C [23 1, which imposes quite stringent constraints for the A-systematics 
in the EKRS analysis, is not used in HKM. As a result, Rp 2 has less shadowing in HKM than in 
EKRS (see also Fig. 1 of 1111 1. especially for heavy nuclei. 

• In addition to the recent DIS data sets, some older ones are used in the HKM analysis. The older 
sets are not included in EKRS. This, however, is presumably not causing any major differences be- 
tween EKRS and HKM, since the older data come with larger errorbars, and are therefore typically 
of less statistical weight in the x 2 analysis. 

• The HKM analysis does not make use of the NMC data [6| on the Q 2 dependence of F 2 Sn /F 2 c . 
As explained above, these data are the main experimental constraint for the nuclear gluons in the 
EKRS analysis. Fig. |3] shows a comparison between the EKRS (solid), the HKM (dotted-dashed) 
results and the NMC data on F^/F^ as a function of Q 2 . As observed there, the HKM results 
do not reproduce the measured Q 2 dependence of F^/F^ at the smallest values of x. This figure 
demonstrates explicitly that the nuclear modifications of gluon distributions at 0.02 < x < 0.1 can 
be pinned down with help of these NMC data. 





(5) 




Fig. 2: The nuclear modifications for the uy, u and gluon distributions in EKS98 (solid) and in HKM (dotted dashed) at 
Q 2 — 2.25 GeV 2 . Upper panels: C/D. Lower panels: Pb/D, with isoscalar Pb. 



3. Future prospects 

3.1 Constraints from pA collisions at LHC, RHIC and SPS 

The hard probes in pA collisions, especially those in LHC, will play a major role in probing the nPDF 
at regions not accessed before. This is sketched in Fig. showing the regions in x and Q 2 probed by 
certain hard processes at different cms-energies. Let us consider the figure in more detail: 

• The measurements of semileptonic decays of D and D mesons in pA will help in pinning down 
the nuclear gluon distributions ll34l l35l . Borrowed from a LO analysis ll35l . the thick solid lines in 
Fig.[2show how the average scale Q 2 = {rrij) of open charm production is correlated with the average 
fractional momentum x = (^2) of the incoming nuclear gluon in dimuon production from correlated 
DD pairs at the LHC (computed for pPb at ^fs = 5500 GeV, 2.5 < y M < 4.0, with no rapidity shifts), at 
RHIC (pAu at = 200 GeV, 1.15 < y M < 2.44) and at the SPS (pPb at </s = 17.3 GeV, < y M < 1). 
For each solid curve, the smallest Q 2 shown corresponds to a dimuon invariant mass M = 1.25 GeV, 
and the the largest Q 2 to M = 9.5 GeV (LHC), 5.8 GeV (RHIC) and 4.75 GeV (SPS) (cf. Fig. 5 in 
[35 1). We observe that data from the SPS, RHIC and LHC will offer constraints in different regions of 
x: the data from NA60 in the SPS will probe the antishadowing of nuclear gluons and the RHIC data the 
beginning of the gluon shadowing, while the LHC data will constrain the nuclear gluons at very small x, 
deep in the shadowing region. 

• Similarly, the DY cross sections of dimuons in pA are expected to set more constraints on the 
nuclear effects of sea quark distributions. In 2 — > 2 kinematics X2 = (M/ ^fs)e~ y , where M and y are the 
invariant mass and rapidity of the lepton pair. The rapidity of the pair is always between the rapidities of 
the leptons forming the pair. The shaded regions in Fig-Qillustrate the regions in x = x^ and Q 2 = M 2 
probed by the DY dimuons in pA collisions at the LHC, RHIC and SPS. The parameters used for the 
LHC are yfs = 5500 GeV, 2.5 y ^ 4.0 with no rapidity shifts [ 36 ]. Note that from the point of view of 
X2 (but not of x\), Pb+Pb collisions at \/ip b p b = 5500 GeV at the LHC are equivalent to pPb collisions 
at yfs pVh = 8800 GeV: the decrease of X2 due to the increase in y/s is compensated by an increase from 

the rapidity shift y |36|: x\ A = (M / Vs pA )e-^ AA+V0 ^ = (M/^ AA )e~ yAA = x£ A . For RHIC, the 
shaded region corresponds to yfs = 200 GeV, and 1.2 < y < 2.2, which is the rapidity acceptance 
region of the forward muon arm in PHENIX [37 ]. For the SPS, we have again used yfs = 17.3 GeV and 
< y < 1. Single muon p^-cuts (see [38 37]) have not been applied, the shaded regions all correspond 
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Fig. 3: The ratio F^/F^ as a function of Q 2 at different fixed values of x. The data is from NMC |6| and the curves are 
EKS98 O (solid) and HKM (T9) (dotted dashed). 

to M 2 > 1 GeV 2 . At the highest scales shown, the dimuon regions at the SPS and RHIC are limited 
by the end of the phase space. The SPS data may shed more light on the question of an EMC effect 
for the sea quarks [39|. The DY dimuons at the LHC will probe the sea quark distributions deep in the 
shadowing region and also at high scales, especially also at Q 2 = M§. The region of RHIC data is again 
conveniently in between the SPS and LHC, mainly probing the beginning of the shadowing region of 
sea quarks. In Ap collisions at the LHC, the DY dimuons with y = (y m - m + y m ax)/2 probe the regions 
plotted with the dashed line (the spread due to the y acceptance is not shown). 

• The dotted line shows the kinematical region probed by open heavy quark production in pA 
(and also in AA), when both heavy quarks are at y = in the ALICE detector frame. In this case again, 
assuming 2 — > 2 kinematics, 

2 = 2my I \f~s aa- To illustrate the effect of moving towards 
forward rapidities, also the case yg = y-Q = 3 is shown. Only the region at Q 2 > m 2 (above the upper 
horizontal line) is probed by bb production. Due to the same kinematics as in open QQ production, the 
dotted lines also correspond to the regions probed by direct photon production. 

• As demonstrated by Fig.[0 the hard probes in pA collisions at RHIC and at the LHC in particular 
will provide us with very important constraints for the nPDF at scales where the DGLAP evolution can 
be expected to be applicable. Power corrections in the evolution BTI l42l can, however, be expected to 
play an increasingly important role towards small values of x and Q 2 . The effects of the first of such 
corrections, the GLRMQ terms I4T1 l42l leading to nonlinear DGLAP+GLRMQ evolution equations, 
have been recently studied in [43 1 in light of the HERA data. From the point of view of the DGLAP 
evolution of the (n)PDF, saturation of gluons takes place when their evolution becomes dominated by the 
power corrections [40 1. Fig-^shows the saturation limits obtained for the free proton and for Pb in the 
DGLAP+GLRMQ approach [44 1 (solid curves; the dotted extrapolations are merely for guiding the eye). 
The saturation limit obtained for the free proton in the DGLAP+GLRMQ analysis should be taken as an 
upper limit in Q 2 but it is constrained quite well by the HERA data (see [43 1). In obtaining the saturation 
limit shown for the Pb nucleus in B4l . the constraints from the Q 2 dependence of F^/F^, however, 
have not yet been taken into account. For a comparison of the saturation limits obtained in other models, 
see l44l . 

• Fig. ^ shows also which hard probes in pA collisions can be expected to probe the nPDF 
in the gluon saturation region. Especially the probes directly sensitive to the gluon distributions are 




interesting from this point of view. Such probes would be open cc production at small px, and direct 
photon production at p? ~ few GeV, both at as forward rapidities as possible (see the dotted lines). Note 
that open bb production is already in the applicability region of the linear DGLAP evolution. In light 
of the saturation limits shown, the chances for measuring the effects of nonlinearities in the evolution 
through open cc production in pA at RHIC would seem marginal. At the LHC, however, measuring 
saturation effects in the nuclear gluon distributions through open cc in pA could be possible. 

3.2 Improvements of the DGLAP analyses 

On the practical side, the global DGLAP fit analyses of the nPDF discussed above can be improved in 
obvious ways. The EKRS analysis should be made more automatic and a proper statistical treatment, 
such as in HKM, should be added. The automatization alone is, however, not expected to change the 
nuclear modifications of the PDF significantly from EKS98 but more quantitative estimates of the uncer- 
tainties and of their propagation would be obtained. This work is in progress. As also discussed above, 
more data contraints should be added to the HKM analysis. More generally, all presently available data 
from hard processes in DIS and pA collisions have not yet been exhausted: for instance, the recent DIS 
data for F% Fe and F% Fe from i/Fe and PFe collisions by CCFRJ451 (not used in EKRS or in HKM), could 
help in pinning down the valence quark modifications [46 1. In the future, the hard probes in pA collisions 
at the LHC, RHIC and the SPS will offer very important constraints for the nPDF, especially for gluons 
and sea quarks. Eventually, the nPDF DGLAP analyses should be extended to NLO perturbative QCD. 
As discussed above, the effects of power corrections I4T1 1421 to the DGLAP equations and also to the 
cross sections [47] should be analysed in detail in context with the global fits to the nuclear data. 
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